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Abstract 



We study the contributions of resonances to D° — D° mixing. Both QQ and 
! hybrid QQG states are considered. Assuming reasonable values for the reso- 

nance parameters, we find relatively sizeable individual contributions to both 
Amo and AT p. We derive a variant of the GIM cancellation mechanism for 



the resonance amplitudes and show that broken SU (3) can allow for appre- 



ciable residual effects. Additional input from meson spectroscopy and lattice 
gauge simulations will be needed to improve the accuracy of these predictions. 

I. BASIC APPROACH 



In this paper, we explore the possible enhancement in the mixing rate of neutral D mesons 
due to nearby resonances. It is interesting that the resonance mechanism, if operative, can 
be available only to D mesons. The light kaon lies below the resonance region and the heavy 
Bd s mesons lie above it. The dynamical mechanism of resonant enhancement constitutes an 
explicit violation of the quark-hadron duality assumption and could influence power counting 
rules built into the HQET estimate of Am D [|],||, which assumes a large energy gap between 
m c and the scale Aqcd at which hadron dynamics is active. In addition, if a resonance is 
viewed as a single-particle intermediate state, then its contribution will be favored over that 
of multibody intermediate states by the 1/N C counting rules. 

For the remainder of this section, we continue the discussion of resonance contributions 
to mixing amplitudes. Section II concerns applications to D° — D° mixing, Section III 
addresses the issue of SU(3) multiplet structure and conclusions are presented in Section IV. 

From standard perturbation theory, the ij th element of the D° — D° mass matrix can be 
represented as 
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The first term in the mass matrix expansion of Eq. (|]) corresponds to the contribution of 
local AC = 2 box and dipenguin operators. These are small in the Standard Model || || . 
Next come the bilocal contributions which are induced by the insertion of two AC = 1 
operators. This class of terms might be enhanced by various nonperturbative effects, and 
therefore is of considerable interest. As follows from Eq. (|lj), one introduces a sum over 
all possible n-particle intermediate states allowed by the corresponding quantum numbers. 
For these continuum contributions, the summation in the second term of Eq. ([]]) takes 
the form of an integral over the energy variable. There will be a unitarity cut in the 
complex energy plane lying along the real axis and beginning at the two-pion threshold. 
The contribution from charged pseudoscalar two-body intermediate states was originally 
considered in Refs. @|| and estimated to be potentially large. However, it remains very 
difficult to reliably determine the total effect associated with n > 2 intermediate states due 
to the many decay modes present, each having unknown final state interaction (FSI) phases. 
For a recent attempt in this direction, see Ref. [0. 

There are also the 'single-particle' effects arising both from bound states and from reso- 
nance intermediate states. A bound state contribution occurs as a pole on the real-E axis. 
A resonance contribution lies in the continuum and corresponds to a pole on an unphysical 
Riemann sheet. Its contribution will be a sharply peaked lorentzian profile of the disconti- 
nuity across the unitarity cut, much like that of a bound state. However, a resonance will 
contribute both to Amo and Ar^. In principle, single-particle effects are rather simpler to 
analyze. The number of such intermediate states is constrained, and they can generally be 
estimated, at least roughly. As mentioned above the l/A c -counting rules, shown already to 
work reasonably well for the estimates of D-meson decay widths, favor a set of single-particle 
intermediate states, i.e. pole diagrams (of Fig. 1) for self-energies of Dl. The light-meson 
and 5-flavored meson single-particle contributions to D° — D° mixing have already been 
analyzed ||. In this paper we study the resonance sector. 

Consider the special role of resonances. To begin, we express the collection of such 
contributions to Am^ (upon neglecting CP- violation) as 
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The pseudoscalar 0~ + (scalar ++ ) intermediate states have CP = — 1 (CP = +1) and 
contribute to the {D$) part of the above. If the mass of the resonance is not too 
far from the D-meson mass, an interesting effect occurs. To highlight the dependence on 
the resonance mass, temporarily consider just the energy denominator in Eq. (|2|). The 
contribution to the energy denominator of a light bound state {e.g. pions or kaons) of mass 
m is P{m 2 ) = l/{m 2 D — m 2 ) = m^ 2 + 0(m 2 /m 2 D ), which amounts to a suppression factor 
of order m^> 2 '. By contrast, the energy denominator (c/Eq. (fj)) for a resonance of mass m R 
and width T R will yield 
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In the limit of a narrow resonance width, the expression for Ar^ becomes proportional to the 
delta function 5{rriR — mo), as expected. The effect of a finite width is to allow the resonance 



2 



to contribute at values itir ^ tub- The contribution to Amp vanishes at m 2 R = m 2 D since it 
undergoes a change of sign there. Considered as a function of the resonance mass thr, the 
maximum effect occurs for m R = m 2 D ±Y R m£> at which Tie P(m 2 D ±T R m£,) = ^fl/(2T R mr>)- 
On the other hand, the resonance contribution to Ar^ is maximized at the different value 
m 2 R = rrijj, but is still substantial at the values m 2 R = m 2 D ± Thw,d which maximize the 
contribution to Amp. In particular, for both Amp and Ar^, the m^ 2 dependence which 
would appear for a very light resonance has been replaced by r^m^. Thus, for a resonance 
sufficiently near the D meson the possibility exists for an enhancement factor of order 
m D/^R — 5 — > 15 for both Amp and AI^ relative to an unenhanced pole contribution. 
Actually it even makes sense to broaden the term 'nearby resonance' to include m R > 1 GeV 
since the m^ 2 suppression mentioned above will be largely overcome. 

Also present in Eq. (fj) are the D-to-R transition amplitudes. To obtain a quantitative 
description, we shall adopt as our AC = 1 hamiltonian the phenomenological Bauer-Stech- 
Wirbel effective operator [ITH], 
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where the constant a-i = —0.55 ± 0.1 is fixed from fitting nonleptonic D decays and it is 
understood that the operator of Eq. (f|) is to be evaluated in vacuum saturation. In vacuum 
saturation, the contribution of resonance R to mixing will be proportional to the squared 
decay constant f R . This has two important consequences: 

1. QQ resonances having J p = + will not contribute, as they occur in P-waves and thus 
have vanishing wave function at the origin. Although they could well be nonzero in a 
more general setting, their absence here suggests they would be suppressed. 

2. We interpret J p = 0~ QQ resonances with masses nearest the D as second radial 
excitations. We include also first radial excitations in our study due to their larger 
decay constants but omit radial excitations above the second. 



II. CONTRIBUTIONS OF INDIVIDUAL RESONANCES 

The fact that the most recent Particle Data Group compilation cites many nonstrange 
and strange resonances in the mass region up to 2100 MeV strongly supports the premise 
of a resonance mechanism. Since experimental data are still relatively sparse in the mass 
region of the D system, however, additional spectroscopic knowledge of this energy range is 
needed. 

In the following, we shall consider the effects of individual QQG or QQ composites. 
Although individual contributions like these will, at least to some extent, be subject to GIM 
cancellations as other states are added in, such contributions nonetheless serve as useful 
indicators of what mixing signal to be reasonably expected. We employ Eq. (|J) for the QQ 
examples, while employing a largely phenomenological method for the QQG case. 

QQ Resonance 

The mixing amplitudes induced by resonance R are 
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where C = 2m D (G F^foCd/ 'a/2) 2 the dimensionless quantities = m 2 R /m 2 D and 7^ = 
Tij/m^i are the reduced squared-mass and width of the resonance, and from the unitarity 
of CKM matrix and neglecting £& we have used £ s ~ — ^ where = V^Vui (i = d,s). 

Since decay widths T ~ 0.20 GeV are characteristic of resonances in the 1 — >• 2 GeV 
mass range, one expects that 1 ^> 7 in applications of Eq. (0). We note that the two mass 
values m res (GeV) = 1.772, 1.973 maximize the resonance contribution to D° — D° mixing 
for fixed decay width and decay constant values. As to the dependence in Eq. (|5]) on the 
resonance decay constant one recalls f n ~ 0.13 GeV, fx — 0.16 GeV for the noncharm 
ground state and fjj ~ 0.2 GeV, fo s — 0.28 GeV PJT2]] for the charm ground state. Ex- 
citations of the constituent quarks will reduce the wave function at the origin and thus 
decrease the corresponding decay constant. For QQ radial excitations, we use hydrogen 
atom wave functions to provide a rough guide in estimating default values of fn. For first 
radial excitations, we estimate fn ~ 0.025 GeV whereas for second radial excitations, we 
use fn ~ 0.01 GeV in our numerical work. We expect contributions from even higher radial 
excitations to be negligible. In our numerical work results are scaled with the square of the 
associated decay constant to allow for any future departures in assumed values or quantum 
number assignments of individual states. 

We display in the Table various resonance contributions (assuming the QQ description) 
to ArriD and to Ar^. The mass and decay width values are taken from the PDG listing. We 
do not intend our listing to be complete, but instead to indicate the magnitudes associated 



with contributions of this type. Although smaller than the current experimental limit p.|TT 
|Am£)| cxpt < 1.3 ■ 10~ 14 GeV, the values are larger than the contribution Amz)| gnd statc ~ 
3 ■ 10~ 17 GeV from the set of pseudoscalar ground state mesons (tt, K, rj, rj'). They also 
tend to dominate the short distance 'box' contributions |Am£>| box — 1.9 x 10~ 17 GeV and 
|Ar D | box ~ 0.75 x 10~ 17 GeV, where we have taken m c = 1.3 GeV, m s = 0.2 GeV, and refer 
effects of QCD radiative corrections to Ref. 0. 

QQG Resonance 

In vacuum saturation, contributions to Am^ from QQ intermediate states arise mainly 
from annihilation amplitudes. It is well known that such amplitudes are subject to helicity 
suppression. This is the same effect which influences the observed patterns of leptonic pion 
and kaon decay modes. It was noted long ago that the effects of helicity supression can 
be (partially) lifted by soft-gluon emission. While this mechanism can be readily applied 
to inclusive heavy meson decay, it is difficult to see how to implement it to the mixing 
matrix elements of Eq. (0) if the intermediate states are of the QQ variety. However, this 
difficulty is avoided if the intermediate state meson is a QQG hybrid state which involves 
a constituent gluon. The argument can be further extended to include penguin operators, 
thus introducing a long-distance counterpart of the dipenguin operator contribution P,|T3 



It is plausible to assume that the gluon produced by a penguin (or any other operator) can 
form a quasibound state along with the uu quark- ant iquark pair, thus producing a hybrid 
resonance (c/Fig. 2). These one-particle intermediate state contributions to Am^ might be 
of importance because of the proximity of the anticipated hybrid-meson mass with 

that of the Z)-meson. 
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In this regard, a particularly interesting candidate is the 7r#(1800) which has J PC = 
h . On the basis of reports from several experimental groups |16| along with various 
quark model analyses [|I7]], it is tempting to assign this particle as a hybrid. One can 
then estimate the contribution of 7r H (1800) to Eq. (§) provided that the mixing amplitude 
g = (D L\H w \ni{(1800)) is known. This amplitude can be inferred from quark models or even 
better, determined phenomenologically by using available data on D decay rates. The idea 
is to search for common decay channels of D and 7Th(1800) where the 7Th(1800) contribution 
is manifest and then estimate the mixing amplitude from this. The situation is as depicted 
in Fig. 3, in which D — 7Th(1800) mixing is followed by 7Th(1800) decay. 



It was noted in theoretical calculations |T7| and hinted at experimentally that the decay 
rates 7Th(1800) — > 7r/o(980), 7r/ (1300) are large for a hybrid 7Th(1800). Thus one can put 
an upper bound on the mixing amplitude g by introducing a model for the resonant decay 
of D-meson via tt h (1800) HT 
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The partial decay width for 7r# — > 7r/ (980) can be written as 

1 2 
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A similar formula exists for the D — > 7r/ (980) transition. In the simplest case, the total 
decay width of tt h can be saturated by the single partial decay width r 7r£r _ +Tr j ( 980 ). This 
is a reasonable approximation as this decay mode becomes dominant for the hybrid tih- 
Thus, using Eq. (0) along with an expression for r £ )^ 7r j ( 980 ), the mixing amplitude g can 
be estimated from 
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Computing the mixing amplitude g using experimental data on the decay rate TE>^nf (980) 
and inserting it into Eq. (0), we estimate | Am^l 7 ^ 1800 ) < 0.3 x 10 -16 GeV, comparable to 
the short distance result. 



III. EFFECT OF MULTIPLET STRUCTURE 

Resonances contributing as intermediate states to D° — D° mixing are expected to occur 
as SU(3) flavor multiplets. The contribution of an entire multiplet will vanish in the limit 
of degenerate light-quark masses due to GIM cancellation. It is not clear how powerful the 
GIM suppression will be, as SU(3) is known to be badly broken in at least some D decays. 
For example, there is the experimentally measured ratio PJl9[] T D o_, K + K - /T D o_, 7T + 7T ~ ~ 3, 



which is unity in the SU(3) limit. Theoretically, it has been suggested that such large 
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breaking is an accumulation of a number relatively minor effects whose ultimate impact is 
substantial. [2(| At any rate, large SU(3) breaking could produce a loophole for evading 
GIM suppression. 

Consider an octet of excited mesons 7r H , Kn, r)n. We use the subscript 'H' to 
represent heavy mesons, and denote the individual members of a resonance octet with the 
above flavor labels. We anticipate the presence of a ninth heavy meson r]' n to allow for 
mixing occurring with rjn. In principle, the mixing angle #h ccn be inferred from either from 
mass determinations or from two-photon branching ratios [j9|,plf. 

We write for the contribution of a mixed octet of resonances to Amp and Tp, 



(10) 

(11) 



where Am^W and AT^ 1 ' are as in Eq. (|5|). The effect of SU(3) breaking can further be 
studied, say for Am^, by expressing the octet decay constant, mass and decay- width factors 
in Eq. (|T0|) as 



fi = fo + Sfi , ^ = Ho + 5 ^ , 7i = 70 + 5% (i = 1, . . . 8) , (12) 

where fo,Ho,Jo an d Sfi, 5Hi, #7i represent respectively the SU(3)-invariant and SU(3)- 
breaking components. This allows for the possibility that the result will be generally influ- 
enced by SU(3)-breaking in the decay constant, mass and decay-width sectors. An expression 
valid to first order in symmetry breaking is 



Am D 
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+^(/o, Ho, 7o) (4<W H - 5 Him ~ 3$Hm) + F -y(fo, A*o, 7o) (^Kh - 5^ - 3Sj m ) 



Ffifo, Ho, To) - S f*n ~ 35 f m ) 
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where for simplicity t/h-^h mixing is ignored and we define 



= 2/i (l -/j ) = a (1 - ho) 2 + (1 ~ 2/i )7g = 2 27q/j (1 - Ho) 

/ " / °(l-Mo) 2 + 7o 2 ' W ° [(l-/,o) 2 + 7of ' 7 " /o [(l-^o) 2 + 7o 2 ] 2 - 

(14) 

Due to the relative lack of data, it is not possible at this time to provide a unique analysis 
of the above relations. Either a small or large effect could emerge, for example: 

1. In the absence of rjn-rj^ mixing, the combination 45hk h - <fyv H - 3<5/%j vanishes by 
virtue of the Gell Mann-Okubo formula and the remaining dependence in Eq. ( p~3"D 
vanishes with the choice Ho — 1 (implying Ff = F 1 = 0). 

2. The choice Ho = 1 — 7o (with 1 ^> 70) yields 

£ f 2 f 2 

J- T 1 J- (/ _ 1 J- " 
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Here the net effect appears in terms of fractional changes in decay constants and decay 
rates, 



{Amp 
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In the first of the above items, the vanishing of Amu| SU( - 3 ^ _brk to first order in SU(3) 
symmetry breaking occurs for a special parameter choice and is clearly more an exception 
than a rule. The second item has the SU(3) degenerate mass set at a lower value (still with no 
r/H-^H mixing) and a nonzero effect will generally occur, although to be more quantitative 
would require additional experimental input. The presence of large symmetry breaking 
might necessitate a treatment beyond the first-order relations given above. 



IV. CONCLUDING REMARKS 

The motivation most often cited in searches for D° — D° mixing lies with the possibility 
of observing a signal from new physics which dominates that from the Standard Model. The 



best experimental limit, recently obtained by E791 [11|, is well beneath existing estimates 



of the Standard Model value. There are plans to improve on the E791 determination, both 



at B-factories [22] and at hadron colliders |23|]. In addition, preliminary plans at Jefferson 
Lab to build a new experimental hall and simultaneously to raise the beam energy suggest 
the possibility for D° — D° mixing studies at that facility as well. For all such efforts, it will 
be crucial to understand the magnitude of the mixing amplitude from the Standard Model. 

In this paper, we have studied the set of potentially significant contributions to D° — D° 
mixing from pseudoscalar resonances. We have shown how an enhancement of order mu /T r 
can arise from a resonance whose mass lies within several decay widths of mo and have also 
have pointed out the importance of lighter resonances due to decay constant dependence in 
the mixing amplitude. In order to obtain a more detailed understanding of the resonance 
scenario, we have considered contributions from both traditional QQ resonances as well 
as an exotic QQG hybrid. Effects of order 10~ 16 GeV are possible for both Amo and 
ATjj. In addition, it would appear possible or even likely in the resonance mechanism that 
\AT D / Arnol — 1 or even larger. This calls into question the usual assumption, that Am D ^> 
ATp, made in searches for CP violation in D-decay using time-dependent measurements []24|| . 
We shall consider generalizations of our approach and implications of our findings vis-a-vis 
CP-violating signals in a separate publication. 

Of course, efforts such as this are severly hampered by a lack of knowledge regarding 
the properties of mesons lying in the 1.6 — > 2.1 GeV mass range. Two kinds of additional 
input would be of significant value to this subject. From experimentalists could come a 
more complete listing of resonance mass and decay-width parameters. Information about 
various decay modes could allow a distinction between the QQ and QQG descriptions. The 
lattice-gauge community could supply numerical estimates of both decay constants of excited 
mesons and also, as a test of vacuum saturation, matrix elements like (clq\Hw\D ) . As such 
valuable information becomes available, we can anticipate real progress in this area. 
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FIGURES 
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FIG. 1. Contribution of resonance R to the D -to-D matrix element. 
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FIG. 2. QQG intermediate state. 
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7T H (1800) 

FIG. 3. Hybrid contribution to D° decay. 
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TABLES 


Table: Magnitudes of Pseudoscalar Resonance Contributions. 


Resonance 


\Am D \ x 10~ 16 (GeV) 


|Ar D | x 10~ 16 (GeV) 


#(1460) 
77(1760) 
vr(1800) 
#(1830) 


~ 1.24 (/x(i460)/0.025) 2 
(0.77 ±0.27) (/,(i 760) /0.01) 2 
(0.13 ±0.06) (/.(i 8 oo)/0.01) 2 
-0.29 (/^ (1830 )/0.01) 2 


~ 0.88 (^ (1460 )/0.025) 2 
(0.43 ±0.53) (A ;(1760 )/0.01) 2 
(0.41 ±0.11) (/.( 1800 )/0.01) 2 
~ 1-86 (/^(i830)/0.01) 2 
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